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This paper studies the electroweak production of all possible four 
fermion states in e + e~ collisions. Since the methods employed to 
\ evaluate the complete matrix elements and phase space are very gen- 

5_i ■ eral, all four fermion final states in which the charged particles are 



detected can be considered. Also all kinds of experimental cuts can 
be imposed. With the help of the constructed event generator a large 
number of illustrative results is obtained, which show the relevance of 
backgrounds to a number of signals. For LEP 200 the W-pair signal 
and its background are discussed, for higher energies also Z-pair and 
single W and Z signals and backgrounds are presented. 
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1 Introduction 



In electron-positron collisions at a few hundred GeV, various gauge-boson 
production processes will be studied. At LEP 200, the reaction 



will be measured, while at somewhat 



W + W~ (1) 

higher energies also 

Z Z , (2) 

Weu e , (3) 

Ze + e~ , (4) 

Z v e v e , (5) 



will come under consideration. The double-gauge boson cross sections de- 
crease, but the single-gauge boson ones slowly increase, with increasing beam 
energy. Around 500 GeV, the cross sections for (H), © and all have about 
the same magnitude, while those for (|2|) and (||) are lower by roughly an or- 
der of magnitude |l| (it should be noted that there exists a number of rarer 
processes in which a single vector boson is produced). 

At LEP 200, the process (HD will be used for an accurate measurement 
of the W mass f2|. Another issue is the test of the non-Abelian couplings in 
reaction (Jj), and, at higher energies, also in the other processes. The effects 
of non-standard couplings between the gauge-bosons have been extensively 
discussed in the literature || . These involve terms additional to the standard 
non-Abelian vertices, as well as altogether new vertices, like ZZj, that would 
also affect processes independent of non-Abelian couplings such as (g) and 
(§• 

It is clear, that for these studies accurate predictions for the cross sections 
(0)"© are needed. Therefore, the effects of radiative corrections should be 
known. Again, a number of such studies is available [|], ||. For W-pair 
production, the radiative corrections have been studied most completely. 
Corrections for the total cross section range from -20 to +3 per cent when 
the collision energy increases from 170 to 500 GeV |3J. 

In addition to the radiative corrections, there is another problem. The 
actually measurable final states in the above reactions are not the gauge 
bosons themselves, but rather their decay products. Thus, all the above 
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reactions are just special cases of 

e + e~ — > 4 fermions . (6) 

To a specific final state, many Feynman diagrams can contribute. Some of 
them are related to the reactions others are not. To distinguish 

them, they will be called signal and background diagrams, respectively. The 
physics issues mentioned above necessitate a good knowledge of the relative 
importance of signal and background diagrams. It may also happen, that dif- 
ferent four-fermion processes lead to the same detectable final state as does 
the signal. Such a process will be called a non-interfering background. 

The aim of this paper is to investigate, at tree level, the contribution of 
signal and background diagrams to the various four-fermion final states. In 
its full generality, this task is complicated for two reasons. 

One reason is that, for a specified final state, the number of Feynman 
diagrams can be quite large. The matrix elements then become complicated, 
and we need an efficient method to compute them. In the present paper, 
spinorial techniques for helicity amplitudes are used. Such techniques 

are most efficient when the fermions are massless. Since we are primarily 
interested in final-state configurations where charged leptons are visible at 
(relatively) large angles to the beams and between each other, collinear situ- 
ations are absent. We therefore assume all the fermions to be massless. Note 
that this implies the absence of diagrams where a Higgs boson couples to the 
fermions; but, although we cannot compute Higgs signals in our approach, 
we can at least reliably estimate the background. 

The second complication arises from the peaking structure. There are 
many different peaks in the multidifferential cross section, each of which finds 
its best description in terms of a characteristic variable. Hence, the number of 
characteristic variables may be much larger than the number of independent 
kinematical variables. We solve this problem by using a multichannel Monte 
Carlo method J3, [KJ. The result of our work, then, is an event generator for 
four-fermion final states. 

The Monte Carlo method allows us to evaluate four-fermion cross sections 
under the imposition of all kinds of cuts. In this paper, we shall do this for a 
number of illustrative examples, indicating the effects of the backgrounds on 
the measured signals. For realistic experimental situations, other effects, like 
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quark fragmentation and hadronization, as well as detector characteristics, 
must be considered. Nevertheless, the examples below already show that 
background diagrams can contribute 4 to 10 per cent to certain four-fermion 
final states. 



It should be noted that the literature contains a number of papers [11 



where also background effects are calculated. None of these considers all 
possible final states. Moreover, our emphasis on obtaining an efficient event 
generator allows for more experimentally precise cross section estimates than 
the results under necessarily simplified phase-space cuts. 



The outline of the paper is as follows. Section 2 lists all four-fermion final 
states for e + e~ collisions, together with the number of Feynman diagrams. 
The next section describes the matrix element computation. In section 4, the 
phase-space distributions are discussed, and illustrative examples of various 
cross sections are given in section 5. In the last section, we summarize our 
conclusions. In an appendix toy model cross sections are given, which serve 
as checks and estimates of phase space integrals. 



2 Signals and Backgrounds 

Since the W ± and Z decay into leptons and quarks, one may classify the 
four-fermion final states into leptonic, semileptonic and hadronic final states. 
The following decays for the W + and Z are considered: 

e + z/ e , , r + z/ T , (7) 

ud , cs , (8) 

e + e~ , , r+r" , v e v e , , v T v T , (9) 

uu , dd , cc , ss , bb . (10) 

In principle, the W knows more quark decays than (|8|). Those that involve a 
b quark are extremely rare, and can safely be neglected. The decay us can, 
usually, experimentally not be distinguished from ud; and the same holds for 
cs and cd. The unitarity of the Cabibbo matrix therefore ensures that the 
use of only (P) will give the experimentally relevant predictions correctly. 
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Figure 1: signal diagrams for e + e~ — * W + W~ . Here, and in the following 
figures, solid arrowed lines stand for fermions, dotted lines denote W's, dashed 
lines Z's, and wavy lines stand for photons. 

The diagrams for the reactions ©-(0) are given in figs. 1-5. In these dia- 
grams, the vector bosons are assumed to be stable. It is convenient to use 
names for certain types of diagrams, as they usually show a characteristic 
peaking behaviour. The first diagram of fig. 1 is called a conversion diagram, 
and the other ones in fig.l are called annihilation diagrams. Thus, ZZ pro- 
duction (fig. 2) consists of conversion diagrams only, whereas Ze + e~ has two 
conversion diagrams. 




\ 

\ 



Figure 2: signal diagrams for e + e — > ZZ. 

In single W production we encounter new types of diagrams. The first and 
second diagram of fig. 3 is called a fusion diagram, the next 4 ones are called 
bremsstrahlung diagrams. The latter represent bremsstrahlung of (in this 
case) a W from Bhabha-like scattering. Similar bremsstrahlung graphs (now 
with a Z radiated off) occur in single Z production, depicted in fig. 4. 




Figure 3: signal diagrams for e + e — ► W + e v e 
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Figure 4: signal diagrams for e + e — * Ze + e 






Figure 5: signal diagrams for e + e — > Zu e iy e 

Note that there is a number of other single W or Z production processes. 
They are of the form 

e + e" -> ^ W + , (11) 

e + e~ -> , (12) 

e + e~ — > fi + fi^ Z , (13) 

e + e" -> z/ M Z , (14) 

e + e~ — > qqZ . (15) 



These reactions receive contributions from the annihilation graphs of figs. 3-5. 
Reactions (|i~3"f ) and (|15|) are also possible through the conversion diagrams of 
fig.4. 
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When the bosons decay, one gets the four-fermion final states. In figs. 1-5 
the decay fermions are attached to the external gauge bosons. In cases with 
identical fermions, more diagrams can arise. For instance, when both Z's 
decay into e + e~, fig.2 will lead to 4 diagrams. Besides the signal diagrams 
of figs. 1-5, there are always a number of background diagrams. The largest 
numbers occur for e + e + e~e~ and e + e~v e v e) the former with only neutral- 
current interactions. The background for e + e~z/ e z/ e contains two bremsstrahl- 
ung graphs, one fusion graph, and 9 multiperipheral diagrams (fig. 6). 




Figure 6: background diagrams for e + e — * e + e v e v e 

For the e + e + e~e~ final state, there are 32 bremsstrahlung diagrams, 4 con- 
version, 16 annihilation, and 32 multiperipheral ones. Of these, fig. 7 depicts 
21 only; the other ones are obtained by permuting the outgoing positrons or 
electrons. 

Since different final states select different sets of diagrams, an inventory 
is made in tables |l|-^| All final states that have the same matrix element are 
grouped together. Moreover, the number of Abelian diagrams (N a ) and of 
non-Abelian diagrams (N n ) is given. We also indicate which of the signals 
(l=(j^),. . . ,5=(|)) lead to the final state, and how many of the background 
diagrams of figs. 6 or 7 occur (iV&). It is seen, that for the leptonic processes 
there are 15 different matrix elements, for the semileptonic processes 10, and 
in the purely hadronic case 7. Since differences in the latter two cases partly 
arise from different quark coupling constants, the numbers of structurally 
different matrix elements in these two cases are 6 and 4, respectively. 
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Figure 7: background diagrams for e + e — * e + e e + e 



3 The matrix elements 

The amplitudes receive contributions from Abelian and non-Abelian graphs, 
with distinct topological structure. They are given in fig. 8. 

In these so-called generic diagrams, all particles are assumed to be outgo- 
ing: assigning two fermion legs to be the initial-state fermions (by crossing), 
the actual Feynman diagrams are generated. In the Abelian diagrams the 
charges of the fermions determine the character of the two exchanged bosons, 
which may be W + , W~ , Z or 7. In the non-Abelian diagrams, two of the 
vector bosons are fixed to be W + and W~ , and the third one can be Z or 7. 
In this way we avoid double-counting of diagrams. 

The particles and antiparticles can each be assigned in six ways to the 
external lines (in principle). This gives, for the Abelian graphs, a maximum 
of 144 different diagrams, and at most 8 for the non-Abelian diagrams. 

We evaluate the matrix element at the level of helicity amplitudes. For 
a specific particle/antiparticle content of the final state, the contributing 
diagrams are enumerated by permuting the assignments of the external legs 
in the generic diagrams, and by relating the charges of the fermions to those 
of the vector bosons V\ and V2. 
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Figure 8: generic diagrams for four-fermion production. The fermion momenta 
and helicities, and the bosons are indicated. The bosons can be either Z, 
W^, or 7; V can be either Z or 7. 

Not all of these generated diagrams will contribute to each helicity am- 
plitude, as can be seen from the following expression for the numerator of 
the abelian diagram of fig. 8: 

= U\{p 1 )^U\{p 2 ) 

X U p (p 3 ) 7^ (^1 + j) 2 + i>Z ) IvUp (p 4 ) 

x u*(p5)Yu a (p 6 ) . (16) 

Here we have disregarded the particle/antiparticle distinction since it is al- 
ready implied by the assignment of the external momenta. The helicity labels 
A, p, a = ± determine the helicity of both external legs on a given fermion 
line. Using the Weyl-van der Waerden formalism for helicity amplitudes || 
(or, equivalently, the Dirac formalism of 0), the expression A can easily be 
calculated. For instance, for A = p = a = 1 one finds 

A(+, +, +; 1, 2, 3, 4, 5, 6) = -4(31>*<46) [(51>*(21> + (53>*(23>] , (17) 
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where the spinorial product is given, in terms of the momenta components, 

by 



(kj) = (p) + ip) 



3 "I 1 / 2 
Pk-Pk 



n ^ 
iPj-Pjj 



(k^j) . (18) 



we denote the expression of Eq . (|17|) by A (l, 2, 3, 4, 5, 6). All helicity combi- 
nations can be expressed in terms of A , as follows: 

A(+ + +) = A (1,2,3,4,5,6) A( ) = 4,(1,2,3,4,5,6)* , 

A{- + +) = 4(2, 1, 3, 4, 5, 6) A{+ - -) = 4(2, 1, 3, 4, 5, 6)* , 

A(+ + -) = 4,(1,2,3,4,6,5) A(--+)=4)(l,2,3,4,6,5)* , 

A(- + -) = 4(2, 1, 3, 4, 6, 5) A{+ - +) = 4(2, 1, 3, 4, 6, 5)* . (19) 

The numerator in the non-Abelian diagrams can also be written in terms of 
the function A: 

u\(pi)j a u\(p 2 ) u p (p-i)-i^Up(p A ) u a (p 5 )j u u a (p 6 ) 
x {g^( Pl + p 2 y + <T (p 5 + p 6 r + 9 U »(P3 + PaT} 
= A(X, p, a; 1, 2, 3, 4, 5, 6) - A(a, p, A; 5, 6, 3, 4, 1, 2) . (20) 

Thus, for massless fermions, every helicity amplitude consists of a sum of 
very systematic, and relatively compact, expressions. When the fermions 
acquire a non-zero mass, the spin states of the external legs are not related 
so simply anymore. Additional helicity amplitudes therefore occur, as well 
as extra terms in the ones above, that are suppressed by one and two powers 
of the ratio fermion mass/energy, respectively. Also, new diagrams crop up, 
where massive fermions are coupled to a Higgs boson. In these diagrams, 
the unsuppressed helicity amplitudes (that nonetheless contain the fermion 
mass in the coupling constant to the Higgs) are precisely the ones that are 
suppressed in the vector boson diagrams, and vice versa. We may therefore 
conclude that the interference of the Higgs diagrams with the vector boson 
diagrams will always be suppressed by a factor proportional to the fermion 
mass squared. 

Finally it should be noted that the vector boson propagators are imple- 
mented in the form q 2 — My + zMyTy, irrespective whether q is timelike or 
not. This only introduces a negligible error. A more refined procedure with 
a q 2 dependent width is of course possible. 
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4 The Monte Carlo 



In this section we describe the structure of EXALIBUR, a Monte Carlo sim- 
ulation program for the four-fermion production processes discussed above. 
First, we shall review the basic ideas, and then show their implementation 
in the actual simulation. 

For the sake of simplicity let us concentrate on the total cross section a 
for some process 

e + (pi)e-(p 2 ) - h(p 3 ) Upa) fs(ps) Up,) ■ (21) 
The main formula is 

a = J f($) rf$ , (22) 

where /(<J>) denotes the matrix element squared (any cut can be easily im- 
plemented by putting /($) = in the unwanted region of the phase space) 
and 

6 

d$ = H d A Pi 5{pf) 5 4 ( Pl +P2-P3-P4-P5- Pa) (23) 

i=3 

is the 8-dimensional massless phase space integration element. 

In order to reduce the variance of the integrand, and therefore the Monte 
Carlo error, it is convenient to introduce an analytically integrable function 
g($), called the local density, that exibits approximately the same peaking 
behaviour of /($) and is unitary, that is, a normalized probability density: 

J g($) d$ = l . (24) 

By multiplying and dividing the integrand by (?($), the cross section 
rewritten as follows 

a — J w($(p)) dp 

where the new integrand 
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is a smoother function of the new set of variables {p{\ defined by 



dp = g($)d$ 
< Pi < 1 (27) 

so that the variance of w(p) is smaller than the variance of /($)• 

When the peaking structure of the matrix element squared is very rich 
(the worst case being e + e~ — > e + e~e + e~ with 144 different Feynman dia- 
grams) one set of new integration variables {pi\ can only describe well a 
limited number of peaks. Therefore a multichannel approach is required in 
which 

N N . 

g($) = Y,a l g l ($) , = ! , 9i ($)d$ = l , (28) 

i=l i=l J 

and where every gj($) describes a particular peaking structure of /($)• Note 
that the conditions on the ctj and <?;($) ensure unitarity of the algorithm, 
i.e. probability is explicitly conserved at each step of the algorithm, without 
additional normalization factors at any stage. 

The numbers ctj are called a-priori weights and, although their numerical 
values are in principle unimportant, they can be used, in practice, to reduce 
the Monte Carlo error [|1(J . 



In EXCALIBUR we have dealt with the problem of the construction of the 
gj($) in a very modular and systematic way. Firstly we have singled out all 
possible kinematical diagrams for every four-fermion final state. Secondly 
we have constructed all building blocks (that is subroutines) necessary for 
the calculation. Finally we have put them together to form the (?«($). 

To illustrate the outlined procedure we shall treat in detail one particular 
channel, that is the conversion channel with one massless and one massive 
particle emitted. 

Referring to the kinematical diagram in figure 9, a suitable choice for the 
8 integration variables is 

• the angle 9 between pi and p^ + p^ in the lab frame; 



^^Here and in the following kinematical diagrams are defined to be pictures that repre- 
sent the various peaking structures. Although they are inspired by the Feynman diagrams, 
they should not interpreted further than that they indicate which variables are most ap- 
propriate to a given pQ |. 
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• the azimuthal angle adjoint to 9; 

• the decay angle 9\ of the particles 3 and 4 in the rest frame of (P3+P4), 
and its adjoint 0i; 

• the decay angle 62 of the particles 5 and 6 in the rest frame of (p$ +pe), 
and its adjoint 02 ; 

• the squared invariant masses s 34 = (p 3 + Pa) 2 and s 56 = (p 5 + p&) 2 - 




Figure 9: kinematical diagram representing a conversion with emission of one 
massless and one massive particle. 

Furthermore we expect in the cross section cos 9 distributed like 

1 , v/i(£ 3 + E A ) - s u 



a — cos ( 



, v ~ 1 , a 



y/s\P3+Pi\ 



S34 like 
and s 56 like 



1 



(29) 
(30) 

\ , (31) 

(s 56 — M 2 ) + M 2 r 2 v ; 

while all other distributions are expected to be more or less flat. Therefore 
the phase space integration can be split into five parts. Using 

l = J ds 3i d A P 3i 5 A (P U -P3- PA>5((p 3 + pi) 2 - s 34 ) , 

1 = J ds 5e d A P 56 <5 4 (P 56 -p 6 - pe) 5((p 5 + p 6 ) 2 - s 56 ) (32) 
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it follows 



l S = I 

J J s 



1/1 "34 

s min ^34 



L 



(v / ^-V / «34) 2 



ds 



56 



[(ss6 - M 2 f + M 2 r 21 



[(s 56 - M 2 ) 2 + M 2 r 2 



X 



x 



n 

(su s 56 \ [ 2w f 1 dcos9 



III 



/•2-7T rl r2n rl 

| J d(f>i J dcosOi | J d(f>2 J dcos9 2 



(33) 



iv v 
where s min is a minimum value for s 34 and 

A(x, y) = 1 + x 2 + y 2 — 2x — 2y — 2xy 



(34) 



The first and the second contributions take care of the photon and the mas- 
sive boson propagators, the third one describes the t-channel distribution of 
cos 9 while the last two integrals represent the two body decays of the mass- 
less and massive boson respectively. The new set of integration variables {p«} 
is 



S34 

s 56 

y + 



cos 6* 

01 

cos 9\ 

02 

cos 9 2 



M 2 + MY tan [p 2 (y + -y~)+y-} , 



tan 



- 1 [(^s-^f-M 2 )/MY_ , 
tan _1 (-M/r) , 

a - [p 4 (a - l)^ + (1 - p 4 )(a + l)^]^ 

2p 6 -l , 
2vrp 7 

2p 8 -l , 



(35) 
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so that it is possible to rewrite 



fd$ = /n dpt sl\ [(«« - M 2 ) 2 + M 2 r 2 ] 
J Jo i=1 ^ ' . ' 

x | A i ^ sfj ( a - cos 



f 7T 7T 

2 2 



0« 9 to 9 V , (36) 

Therefore the total local density can be written as a product of five contri- 
butions 

9($) = 9i ■ 9u 4 9m ■ 9iv ' 9v • (37) 

The advantage of this splitting is that every part of the algorithm can be used 
again for other channels, both to generate distributions through eqs. (|35|) and 
to compute local densities. For example, using the same ingredients, the 
channel in figure 10 can be easily built as follows 

<«* - r [(856 _ M t; +M2n -K^-^ + M^] 



x 



u 

<^456 ui 
Vl ' S 456 
s 56 *456 



x « A2 i7— H d(p L a-cos^ -^-^' 



/■27T /-l 

X 



/•27T ^1 /*27r ^1 

| / d0i / dcosOi | / rf^2 / d cos 9 2 (38) 

i/0 k — 1 " J — 1 



IV 



where s 3 = p\, s 456 = (p4 + P5 + Pe) 2 and (#1, 0i), (9 2 , 02 ) are the dacay 
angles of S456 and S56. 
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Pi 




P2 



Figure 10: kinematical diagram representing a decaying massive boson radiated 

off a final leg. 

This finishes our description of the event generation procedure. More 
details can be found elsewhere ||16|| . Finally, a brief remark is in order on the 



collinear singularity in the multiperipheral and bremsstrahlung diagrams. 
Consider an incoming electron with momentum p M , which scatters over a 
small angle into an outgoing momentum g M , by emitting a photon with vir- 
tuality t. We then have 

t = (p - qf = 2m\ ~ 2 A° + %\p\\q\ cos9 (39) 

where 9 is the scattering angle. For very small scattering angles, the cross 
section will indeed be dominated by the Feymann diagrams with this multi- 
peripheral structure. If the electron is massive, there is no singularity, in the 
sense that 

t < t (40) 

to = 2(ml - A° + \p\\qi) ~ - mliP ° ~ %f (41) 

PoQo 

where we have assumed m e << qo,Po- The total cross section is therefore 
finite. For massless electrons this is not the case because the collinear sin- 
gularity can be reached. We may, however, mimic the dominant effects of a 
nonzero mass by imposing, on the generated event, a cut 

> 9 (p ,q ) (42) 

0o(Po,go) = — — • ( 43 ) 
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This cut leads, for massless electrons, to the same limit, to, on t that would 
arise in the massive case. In our Monte Carlo approach, such a cut can be 
easily applied to one or to both multiperipheral legs. In fact, in the t-channel 
generation of 9, the value of the exponent v in eq. |29] is arbitrary. Therefore, 
using v < 1, cos 9 can be generated in its full range [-1,1] even in the case 
of massless particles and any final state dependent cut implemented in an 
event by event basis. Note, however, that this procedure only reproduces the 
leading log result of the full, massive, cross section. For instance, subleading 
terms due to spin-flip amplitudes are necessarily absent. On the other hand, 
the angle 9q is typically much smaller than any detector is likely to cover, so 
the kinematics of the final state are, for all practical purposes, those of the 
massless electron cross section. 



5 Results 

In this section we present some illustrative examples of the physics results 
one can obtain with the event generator. The input parameters used are 
the masses and widhts of the gauge bosons and the electroweak couplings 
parametrized by a and sin 2 9w- The masses and the widths are taken to be 
independent parameters. The actual values in the program are a = 1/129, 
sin 2 ^ = 0.23, M w = 80.5, T w = 2.3, M z = 91.19 and T z = 2.5 (all 
GeV). Note that, for simplicity, we have obtained results for fixed, energy- 
independent widths. The change to s-dependent widths is trivial by a simple 
modification of the boson propagators computed for each Monte Carlo event. 
There are three types of results. The first type consists of checks against other 
calculations. The second type is concerned with W-pair production at LEP 
200 and somewhat higher energies. The third type is related to high energy 
results for a future linear e + e~ collider and illustrates signal/background 
issues for all five signals. 

As to the checks of the program, we first have verified that all toy model 
matrix elements described in the appendix are, in fact, correctly computed 
by the program. This serves to establish the consistency and reliability of 
our approach. As to more physically meaningful comparisons, results of 
the generator for the W-pair signal have been obtained for a number of 
decreasing widths. The extrapolated values compare well with the zero width 
analytic result (table f|). In order to compare with a more complicated set of 
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diagrams the semileptonic channel /i v^ud was chosen for which an analytic 
calculation exists Both the full calculation with all diagrams and the 



signal calculation agree within the errors of the Monte Carlo calculation. 

The LEP 200 results are divided up into leptonic, semileptonic and hadronic 
final states. 

As for the leptonic final states, results for two charged leptons and two 
neutrinos are collected in table |5|. The first entry is the W-pair signal, then 
the final states 1,6,16,20,18,2 and 8 of table [I] follow. They are calculated for 
all diagrams. Reactions 16,18 and 20 are the non-interfering backgrounds. 

Focussing on the results for a c.m.s. energy of 200 GeV, we see that taking 
into account all diagrams for the final states e + e~v e u ei e~b , e v fl {i + , /i + ji~ v^v^, 
/ u~z/ M z/ r r + gives generally speaking a larger cross section than the W-pair 
diagrams alone. One finds increases of 12, 3, 5, % respectively. The non- 
interfering backgrounds to the e + e~ final state add up to 9.5 % of the WW 
signal and for the [i + or r + r~ final state to 9%. 

For completeness we note that there is also an invisible 4 neutrino cross 
section. Without any cuts it has at 200 GeV the value .8035 10 _1 pb, to be 
compared with the total 2 neutrino cross section of 41.82 pb. 

The semileptonic states related to the W-pair signal are final states 1 
and 5 of table |[ The results involving only the signal diagrams and those 
containing all diagrams are listed in table [| Inclusions of all diagrams for 
the e~v e ud final state gives an increase of 4 % above the WW signal at 200 
GeV. When an invariant mass cut is imposed on the quark pair, the effect 
does not change, while a double invariant mass cut (provided the lepton- 
neutrino invariant mass can be experimantally reconstructed) washes out 
any difference, as shown in table [7|. For the ^Tv^ud final state there is little 
difference between the signal and the full calculation at LEP 200 energies. 

The results for the 4 quark final states are given in table The full 
calculation for the uudd final state gives an answer 3% above the signal value, 
again at 200 GeV. For the udsc there is no difference between the signal 
and the full calculation. When one adds all 4 jet events coming from the 
signal and does the same for all non-interfering backgrounds, then the latter 
are 8% of the real signal. Note that those are non-interfering electroweak 
backgrounds and not QCD backgrounds. These should be added as well, but 
this is beyond the scope of the present paper. 

The third type of results is for energies of 500, 1000 and 2000 GeV. 
It should illustrate the possibilities of extracting the signal from a specific 
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final state, where the final state is produced by all diagrams. We take the 
idealized case where one could measure any invariant mass combination of 
the four fermions in the final state. When a charged or neutral fermion pair 
has an invariant mass m (ij) in the interval [My — 2r, My + 2r] for V= W or Z 
we identify the fermion pair as originating from the vector boson V. In table 
H the e~z/ e z/ e e + final state produced at 500 GeV is analyzed with respect to 
its vector boson content. There are events which do not contain any mass 
combination giving rise to a vector boson, but there are also some events 
containing apparently three vector bosons. It is now interesting to compare 
this analysis carried out on atot containing all diagrams with a similar analysis 
on <j(i), the cross section originating from the signal diagrams for process (i). 
Tables |TU] - do this for 4 different final states. In table [T^ the e~z/ e z/ e e + 



final state is analyzed. The cross section a to t always is greater than the sum 
of all signals. This demonstrates the relevance of the diagrams of fig. 6. 
Tables |11| - [13] consider e~v e ud, e + e~uu and u e u e uu final states. Usually 



a to t is greater than X^c(t), but not always. When one applies the idealized 
search algorithm to a to t and to every a{%) both a to t and a(i) are reduced, 
the former of course much more drastically. The vector boson cross sections 
extracted from atot tend to be greater than those originating from the signal 
cross section with a possible exception at 2000 GeV. The extracted WW 
cross sections are within the errors equal, for Ze + e~ and Zvv this is also 
more or less the case at 1000 and 2000 GeV. 

From these idealized examples it seems that search algorithms are possible 
which extract the specific vector boson signal cross section. This causes a 
significant loss in statistics. Since most of the theoretical anomalous vector 
boson coupling studies have been performed for non decaying vector boson 
final states, one has to take into account this reduction of statistics. When 
one prefers to keep the higher statistics, then it seems recommendable to 
study the effect of non-abelian couplings directly on the four-fermion final 
states. 

It should be noted that the energy dependences of the signals in tables 
[It] - [T3] can differ from those in ref. ||T|. The reason is that in the tables the 
cross sections always have some cuts on the scattering angle, whereas in [|TJ 
the full phase space has been taken into account. 
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6 Conclusions 



In this paper one strategy is chosen for the computation of all possible four 
fermion final states. Whatever process is chosen the calculation proceeds 
through the same series of steps. On one hand, this concerns the matrix 
element which in principle contains all diagrams obtained by permutation 
from two basic diagrams. On the other hand, it means that, for a specified 
final state, the possible peaking structures are found and the phase space is 
generated accordingly. The algorithms for the generation of phase space are 
always made up of a particular combination of elementary building blocks. 
The use of one procedure for the matrix element and of a set of well tested 
(through model matrix elements) algorithms for distributions of certain kine- 
matical variables makes the calculation transparent and less error prone. The 
neglect of fermion masses renders the computation fast when one considers 
the often large number of Feynman diagrams. 

The four fermion final states offer a rich phenomenology: quite different 
pure leptonic, semi leptonic and quark final states can be studied. Moreover, 
the number of contributing diagrams to a specific process varies a lot. In 
terms of produced vector bosons, which subsequently decay, there are five 
signals which have been treated frequently in the literature. For these the 
questions of signal and background are discussed in this paper for a set of 
four fermion final states. It is shown that sometimes the backgrounds are of 
the same order of magnitude as radiative corrections and that certain exper- 
imental cuts can reduce the backgrounds. Often both signal and background 
are then reduced, which causes a loss in statistics. Finally, some possible 
future extensions of our calculation or their applications are listed 

1) certain studies, like effects of anomalous couplings have actually been 
performed for vector boson final states. In view of the above losses 
in statistics one may consider to build in directly anomalous couplings 
into the matrix element of section 3; 

2) another possibility is to use the event generator for simulating non- 
interfering background events for Higgs searches; 

3) when the QCD matrix elements are also taken into account, a full 
analysis of four jet events in the region of vector boson production can 
be performed; 
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4) radiative corrections to signals and backgrounds have not been con- 
sidered here. Some effects could be included in a further development 
of the event generator. The most important one is the initial state 
radiation, either in one loop or in a leading logarithmic approximation. 

A last remark is in order. When one would like to study the precise effects 
of untagged electrons or positrons one clearly needs massive fermions in the 
matrix element and also an efficient algorithm for generating two photon 
multiperipheral diagrams in the full phase space. This problem has been 
dealt with in the literature in the past. It seems that one is unavoidably led 
to time consuming event generators. Therefore an extension of the present 
event generator to the exact massive fermion calculation is bound to be quite 
involved and conducive to an appreciable loss in program speed. 
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Appendix A: Model matrix elements and cross 
sections 

It is obvious that a Monte Carlo program such as EXCALIBUR has to be tested 
extensively before it can be considered reliable. One of the most important 
kinds of tests is the Monte Carlo integration of a matrix element that, at 
the one hand, displays some of the salient features of the cross section, and, 
on the other hand, has an analytically known cross section. Note that this 
integral can be computed for any allowed set of a-priori weights cej, that 
is, either with the Monte Carlo channel that gives the same peaking 

structure, or with another channel, or with any combination of channels, and 
in each case the correct cross section should be reproduced (with of course 
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a larger error if the appropriate channel is not included). Therefore, we give 
in this appendix a list of such model matrix elements. These can be used 
either for a user test of the Monte Carlo, and also (when the relevant coupling 
constants are added) to get a rough idea of the order of magnitude of the 
contribution of the various kinematical regions to the total cross section. In 
all cases, pi )2 denote the incoming momenta, and ^3,4,5,6 the produced fermion 
momenta. Each model matrix element (squared and spin summed/averaged) 
is denoted by Tj, and the corresponding total cross section is 

A.l Conversion channels 

We distinguish the three cases of conversion of two photons, one photon and 
one resonance and two resonances, respectively. 



Pi 




P3 Pi 




P3 Pi 

Pa 

P5 



P3 

Pa 
Pb 
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Figure 11: Figure 12: Figure 13: 

two-photon photon-resonance two-resonance 
conversion conversion conversion 

For these channels we define the function 

-1 



X(s,s 1 ,s 2 ) 



log 



(< 



s — Si 



s- si - s 2 + A(s, si, s 2 ) 
\(s,si,s 2 ) 

1/2 



si - s 2 

\2 



S2 



Asis 2 



(45) 



This function helps us to obtain analytical closed forms for the various cross 
sections. Owing to the logarithm, it depends not too strongly on si and s 2 . 
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A. 1.1 Two photons 

The kinematical situation is depicted in fig. 11. We have s± — qf — G03+P4) 2 , 
s 2 — Q 2 = (j?5 + Pg) 2 an d = Pi ~ <?i ■ I n order to avoid singularities we 
impose a lower bound on Si and s 2 , and our model for Ti becomes 

T = r(su s ) QJsy > gg) 



and the cross section 

3 



(47) 



Here and in the following, the step-function of an inequality is 1 if the 
inequality holds, and otherwise. Furthermore, w = y/s, Wi = y/sl (i= 
0,1,2) and F = j^)s- 

A. 1.2 One photon and one resonance 

One of the photons now becomes a Breit-Wigner resonance with mass m and 
width T (see fig. 12). We may now drop the lower limit on si. 

ST(S 1 ,S 2 ) 0(S 2 >S ) , . , 2 . 2 

T2 = — \Fa\ ^T~\ ' R{si) = \si-m + imY\ . (48) 

\Q 2 \ i?(si) s 2 

The total cross section is given by 

^V 3 s _ 2 



° 2 = F UJ 2unT [MZ) ~ MZ * )] ' z = m +tmT 
A 2 (z) = 2[B 2 (y) + B 2 (-y)] , y = yfl 

B 2 (y) = log (^) log (^LZl) 
\ w J \ -y ) 

_ u fw^^A t^v\ (49) 

V / \w-yj 

Evaluation of this cross section relies on the evaluation of logarithms and 
dilogarithms of complex arguments. Since their imaginary parts are not 
infinitesimal, this poses no principal difficulty 
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A. 1.3 Two resonances 

We use, for generality, two different masses m 12 and widths r 12 . Our model 
matrix element squared is now 

T 3 (s 1: s 2 ) = p , v p / v , ^fe(sfe) = |sfe-m^ + im fc r fc | 2 (50) 

/tl(SiJK 2 (S2j 

and the total cross section a 3 is quite analogous to er 2 ; 

7T\ 3 — S 2 

0"3 



2/ 4mim 2 rir 2 
x [A 3 (z!,z 2 ) - A 3 ( Zl ,z*) - A 3 (z*,z 2 )+A 3 (z*,z*)} , 

z k = m 2 k + im k T k , 

A 3 (z 1: z 2 ) = B 3 (y 1: y 2 ) + B 3 (y 1: -y 2 ) + B 3 (-y 1: y 2 ) + B 3 (-y 1: -y 2 ) , 



X 



log I 1 - — I + log I — I - log '"' " ! 



yij \yi + y2-wj \yi + y2-w / 

logfl-^logf^-logfl-^lo^ V ' 



yij \wj v w) \yi + y2-wj 

+ log (-*») log (-^-) • (51) 
V wj \yi + y2-w) 

A. 2 Non-Abelian boson fusion channels 

Since the fusing bosons are in the t-channel, we neglect their Q 2 dependence 
if they are W or Z. This approximation is not so good if s >> m 2 zw , but 
at LEP 200 energies it is quite accurate and, in any case, we are allowed to 
propose any form for T, provided it can be checked with our Monte Carlo. 



23 




Figure 14: two Figure 15: two Figure 16: boson 

heavy bosons fusing heavy bosons fusing and photon fusing 
into a photon into a resonance into a resonance 



A. 2.1 Fusion of two heavy bosons into a photon 

The kinematics is given in fig. 14 and we write si = (p^ + Pi) 2 - Our model 
for this channel is 

- - ^ > So) (52) 



and the cross section evaluates to 

3 



2£o\ ^ jS 5 2sq 

s J s 2 s 2s 2 



(53) 



A. 2. 2 Two heavy bosons fusing into a resonance 

Replacing the photon of fig. 14 by a resonance, we get fig. 15 and we have 

1 



and 



Aziz) 



R(si) 



, R(si) = \s\ — m + imT\ 



(54) 



,27 AimT 



[A 5 (z) - A 5 (z*)] , z = m 2 + imT 



5 
2 



z 2z T /s 
- + — Li 2 - 

s s \z 



(55) 



A. 2. 3 One photon and one heavy boson fusing into a resonance 

To avoid collinear singularities, we now have to impose an angular cut. In 
fig. 16, let ip denote the angle between pi and p*5 in the lab frame. We impose 
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a symmetrical cut < ip < it — ipQ. Modelling the (complicated) photon 
propagator by a simple angular dependence, we may put 



6{ip < ip < vr - Vo) 



1 — cos tp 

So the cross section is obtained immediately: 



1 (1 + cosV>o \ ~ 

0-6 = 77 log " CT 5 

2 V 1 — cos tpo J 



(56) 



(57) 



A. 3 Annihilation channels 

Since the total energy is fixed, it is immaterial whether the boson in which 
the e + e~ annihilate is a photon or a resonance (of course, for the above 
mentioned order of magnitude estimates one may want to put in the correct 
energy dependence). 
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Figure 17: annihilation Figure 18: annihilation 
into two resonances into a photon and a 

resonance 




Figure 19: annihilation Figure 20: annihilation 
into a photon and two into a resonance and two 
fermions fermions 
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A. 3.1 Annihilation into two resonances 

The kinematics is given in fig. 17 and, as before, Sj = q 2 . We now take 



2 

S 



r, = *(«,,«,) fli(8i)fl2(S2) ■ (58) 

where we have introduced a factor analogous to r(s 1; s 2 ): 

X(«i> s 2) = T7 r 0(tt>i + w 2 < ^fs^) , s m < s . (59) 

A(s, si, s 2 ) 

The cutoff is intended to avoid the singularity in A at w\ + w 2 = w. The 
integration is very similar to that of T 3 : 



^ = F - 



7T\ 3 — S 2 



,2/ 4mim 2 TiT 2 
x [A 7 ( Zl , z 2 ) - A 7 ( Zl , 4) - A 7 (z*, z 2 ) + A 7 {z\, 4)] , 

A 7 (^,z 2 ) = [A 3 (^ 2 )U Sm . (60) 

A. 3. 2 Annihilation into a photon and a resonance 

This channel can, strictly speaking, not occur in e + e~ collisions, but we 
include it for completeness, with an eye to, for instance, electron-quark scat- 
tering. Fig. 18 gives the kinematics. 
We put 

s fl(si > 0) 

T8 = X(S1 ' S2) Sl R(s 2 ) (61) 



7T \ 3 S 



so that we have 

A B (z) = [A 2 (z)],^ . (62) 

A. 3. 3 Annihilation into two fermions and a photon 

This situation is given in fig. 19. A cut on si = {pz + P4) 2 = q\ now 
automatically imposes a cut on s 2 = (p$ + q\) 2 , and we may use 

0{s l > s ) (63) 



SlS 2 
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leading to 



09 



2 + £°)logi- + 3(l-^ 
s J s V s 



(64) 



A. 3. 4 Annihilation into two fermions and a resonance 

Replacing the photon of fig. 19 by a resonance, we obtain fig. 20. We have 
the following model: 

Tie 



and 



^io0) 



7T 



F 

3-2(1 



2 7 2imT 

z 



s 2 R(si) 

[A w (z)-A 10 (z 
log (l - 



1 



f ) Li = (f 



A. 4 Final state bremsstrahlung channels 



(65) 



(66) 



In all these boson is exchanged in the t channel and a fermion pair 

radiated off an outgoing particle. 




Figure 21: photon Figure 22: resonance 

bremsstrahlung in heavy bremsstrahlung in heavy 
boson exchange boson exchange 
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Figure 23: photon 
bremsstrahlung in 
photon exchange 



Figure 24: resonance 
bremsstrahlung in 
photon exchange 



A. 4.1 Heavy boson exchange 

If the radiated fermion pair comes from a photon, the situation is given in 
fig. 21. Since there is (by the definition of our model) no angular dependence 
arising from the boson exchange, we simply have 

Tn = T 9 , on = cr 9 . (67) 

On the other hand, for resonance bremsstrahlung, the appropriate model for 
fig. 22 is again one that we have already mentioned: 

Tu = Tio , 5-12 = 5"io • (68) 



A. 4. 2 Photon exchange 

Now, we have a nontrivial angular dependence, which we again simply model 
by the scattering angle ip between p% and ps, in the lab frame. For photon 
bremsstrahlung (fig. 23) we have a simple modification of Tn; 

6{^ Q < i> < vr - ipo) 



T 



13 




1 — cos ip 



. . 1 + cos^ \ _ ( . 

(Ji3 = -log cr u . (69) 



1 — COS ljj J 

Similarly, for the resonance bremsstrahlung of fig. 24 

9(i/j < tp < n - 4> 



T 



14 




1 — COS 1p 

f cos ip 

1 — COS -00 



. . 1 + cos^ \ _ , . 

au = - log ai2 . (70) 
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A. 5 Initial state bremsstrahlung 

The fermion pair is radiated off an initial leg. This complicates matters 
somewhat. 




Figure 25: photon Figure 26: resonance 
bremsstrahlung in heavy bremsstrahlung in heavy 
boson exchange boson exchange 




Figure 27: photon 
bremsstrahlung in 
photon exchange 



Figure 28: resonance 
bremsstrahlung in 
photon exchange 



A. 5.1 Heavy boson exchange 

Since the heavy boson propagator is again assumed to give a constant angular 
distribution, these channels look somewhat like the conversion channels. We 
define, for fig. 25 (si = q\ = (p 3 + p A ) 2 , s 2 = (p 5 + p 6 ) 2 , Q M = Pi - <?i ) 



which leads to 



Cl5 



15 



i r 



t(si,s 2 ) 6>(si > s ) 



log hi 

s s 



4 1 



Wo 

w 



(71) 



(72) 
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For fig. 26, we have, in analogy: 

m t(si,s 2 ) 



16 



\Q 2 \R(si) 



(73) 



and 



A 16 (z) = B 16 (y) + B 16 (-y) 

B 16 (y) = ( w - y )Hog(l--) - — + wy . (74) 

V y J y 

A. 5. 2 Photon exchange 

For these channels, the situation is more complicated. It is straightforward 
to define the scattering angle ip between p 2 and p Q not in the lab frame, 
but in the rest frame of ps + p$, but a simple cut on this angle translates 
to a very complicated one in the lab frame. We therefore again rely on 
our freedom to postulate any reasonable behaviour, and simply model the 
photon propagator such that it is simple in terms of this new if), and non 
singular upon integration. At small scattering angles, this is actually not 
a bad approximation to the real situation. For photon bremsstrahlung (fig. 
27) we therefore propose 

T17 = T 1B j 1 —— , a < 1 . (75) 

(1 — cos ip) 



So that 



^=[T^r i5 - (76) 



Similarly, for the resonance bremsstrahlung of fig. 28, we use 

1 



Ti 



18 — 1 16 



(l-cos^) a ' 



*18 = ( j^J -16 • (77) 

Note that the algorithms used in the Monte Carlo allow us to actually gener- 
arate such non-singular angular distributions without restriction on the scat- 
tering angles. 
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A. 6 Multiperipheral scattering 

We can distinguish according to the type of exchanged bosons 
Pi v ^ P3 Pi \ ^ P3 Pi 



P2 




P2 




Figure 29: Two heavy Figure 30: photon and Figure 31: Two 
bosons fusing into heavy boson fusing photons fusing into 
fermions into fermions fermions 

We shall neglect the boson propagators for simplicity. Again we are free to 



do this. We have s\ 



(P3 + Pa) , s 2 = {ps+Pe) , 



Pi 



the exchange of two heavy bosons we have simply the model 
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t(si, s 2 ) 
s 2 \Q\ 2 



019 



F{ 1YL 

.2 J 6 



gf. For 



(78) 



Now for each photon, we apply the trick mentioned above. Let if>i be the 
angle (pi,P3) in the (P3+P4) rest frame and ip 2 that of {p 2l po) i n the (ps+Ps) 
rest frame. Then, for fig. 30 and 31 respectively, we write 



T 2 q 
T 21 



Ti 9 
T 2 q 



1 



^1 - COS^l)" 

1 



020 



021 



(1 - cos ip 2 ) a 
This finishes our list of matrix element models. 




019 



020 



(79) 
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Table 1: leptonic four-fermion final states in e + e collisions. 
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label 


final state 


N a 


N n 


total 


N b 
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e~v e ud 
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20 
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48 


26 


2 4 


14 


e + e~cc 
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e + e~dd 
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48 


26 


2 4 
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e + e~ss 












17 


e + e~bb 
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fI + fi~UU 


24 





24 


10 


2 


19 
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T + T~UU 












21 


T + T~CC 













Table 2: semileptonic four-fermion final states in e + e collisions (continued 
on next page). 
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label 


final state 


N a 


N n 


total 


N b 


signals 
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Table 2: semileptonic four-fermion final states in e + e collisions (continued 
from previous page). 
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label 


final state 


N a 




total 


N b 


signals 


1 


uudd 


33 


2 


35 


2 


1 2 


2 


ccdd 












3 


udsc 


9 


2 


11 





1 


4 


dues 












5 


uuuu 


48 





48 


20 


2 


6 


cccc 












7 

8 


dddd 
ssss 


48 





48 


20 


2 


9 


bbbb 












10 


uuee 


24 





24 


10 


2 


11 


uuss 


24 





24 


10 


2 


12 


uubb 












13 


ccdd 












14 


cebb 












15 


ddss 


24 





24 


10 


2 


16 


ddbb 












17 


ssbb 













Table 3: hadronic four-fermion final states in e + e collisions. 





yfi = 165 GeV 


yfi = 175 GeV 


v /i= 185 GeV 


^ = 195 GeV 


An. 


11.95 


18.48 


20.41 


20.82 


M. C. 


11.96±0.03 


18.45±0.04 


20.40±0.05 


20.79±0.05 



Table 4: analytic vs. Monte Carlo W-pair signal (in picobarns). 
Here a = 1/129, M w = 80.5 GeV, M z = 91.9 GeV, cos9 w = M w /M z and 
T w = (3aM M /)/(4sin 2 Q w ). The zero width extrapolated values of M.C. are 
obtained with the substitutions T\y — > T^/N and a — > u/N 2 , by taking the 
numerical limit for growing N. 
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process 


^=150 GeV 


y/s = 175 GeV 


y/s = 200 GeV 


yfe = 500 GeV 


W-pair 


.3600 10- 2 
.0011 


.1181 
.0002 


.1304 
.0003 


.2130 10- 1 
.0011 


e + e~u e u e 


.2949 10~ 2 
.0008 


.1208 
.0002 


.1466 
.0003 


.6111 10" 1 
.0028 




.4350 10~ 2 
.0012 


.1196 
.0002 


.1364 
.0003 


.2209 10" 1 
.0011 




.8379 10- 3 
.0026 


.1906 10~ 2 
.0004 


.6226 10~ 2 
.0008 


.2342 10~ 2 
.0006 


v T u T n + n~ 


.7927 10~ 3 
.0031 


.1342 10" 2 
.0003 


.5457 10~ 2 
.0007 


.1054 10~ 2 
.0003 




.1002 lO" 2 
.0004 


.1842 10~ 2 
.0005 


.6004 10~ 2 
.0009 


.7660 10~ 2 
.0079 


W-pair 


.3598 10~ 2 
.0011 


.1187 
.0002 


.1315 
.0003 


.2130 lO" 1 
.0011 




.2912 10~ 2 
.0008 


.1192 
.0002 


.1359 
.0003 


.3598 10- 1 
.0014 


(x-is fl is T T + 


.3649 10~ 2 
.0010 


.1193 
.0002 


.1307 
.0003 


.2101 lO" 1 
.0010 



Table 5: leptonic cross sections in picobarns. The second line of each entry 
is the Monte Carlo error. In the first six entries m(e+e-),m (/*+/*-) > 10 GeV , 
| cos# e ± ; M ±| < 0.9 (9 is the scattering angle), E e ± t M ± > 20 GeV. In the last 
three ones E e - t (U ± ; T + > 20 GeV, | cos# e - ; M ± t T +| < 0.9. 
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process 


yfs = 150 GeV 


y/s = 175 GeV 


y/s = 200 GeV 


^ = 500 GeV 


W-pair 


.9605 10~ 2 


.3342 


.3601 


.4682 lO" 1 




.0030 


.0007 


.0008 


.0028 


e~ v e ud 


.7563 10~ 2 


.3345 


.3745 


.8293 10- 1 




.0021 


.0007 


.0009 


.0035 




.9728 10~ 2 


.3337 


.3610 


.4599 10- 1 


.0028 


.0007 


.0008 


.0024 



Table 6: semileptonic cross sections (in pb). E e - „- u £ > 20 GeV, 
| cos9 e - > „- Uj j| < 0.9 , | cos l(u,d) | < 0.9, m(ud) > 10 GeV. 





y/s = 190 GeV 


y/s = 200 GeV 


yfi = 500 GeV 


W-pair 


.3532 


.3419 


.2498 lO" 1 




.0009 


.0008 


.0018 


1 cut 


.3612 


.3550 


.3522 lO" 1 




.0009 


.0009 


.0020 


W-pair 


.3392 


.3259 


.2236 lO" 1 




.0009 


.0009 


.0014 


2 cuts 


.3407 


.3292 


.2228 lO" 1 




.0009 


.0009 


.0014 



Table 7: e (3) v e (4) u (5) d (6) cross section in picobarns. The l si entry 
is a comparison between the WW signal and the full result with a cut 
70 < m {56) < 90 GeV. In the 2 nd entry 70 < m {34) , m (56) < 90 GeV. All the 
other cuts are like in table |], except that an energy cut E e - u j > 50 GeV 
is required at = 500 GeV. 
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process 


v /i= 150 GeV 


v /i= 175 GeV 


y/s = 200 GeV 


yfe = 500 GeV 


W-pair 


.2141 10- 1 
.0019 


.7699 
.0018 


.8726 
.0023 


.1083 
.0008 


uudd 


.2317 lO" 1 
.0010 


.7745 
.0018 


.8987 
.0023 


.1103 
.0006 


udsc 


.2132 10- 1 
.0010 


.7720 
.0018 


.8709 
.0023 


.1052 
.0006 


uuuu 


.1033 10~ 2 
.0008 


.1956 10- 2 
.0006 


.1031 10- 1 
.0002 


.1845 10-2 
.0006 


dddd 


.3900 10~ 3 
.0018 


.1171 10-2 
.0003 


.1508 10- 1 
.0003 


.2239 10-2 
.0008 


uucc 


.2417 10~ 2 
.0012 


.4032 10-2 
.0014 


.2066 10- 1 
.0004 


.3715 10-2 
.0013 


uuss 


.1836 lO- 2 
.0009 


.3474 10-2 
.0011 


.2499 10- 1 
.0005 


.4143 10-2 
.0014 


ddss 


.8245 10~ 3 
.0038 


.2437 10-2 
.0007 


.3026 10- 1 
.0006 


.4475 10-2 
.0017 



Table 8: hadronic cross sections (in pb). E( a u particles) > 20 GeV, 
I cos 9 ( a u particles) \ < 0.9. Moreover m(ij) > 10 GeV and |cosZ(«,j)| < 0.9 
between all possible final state couples. 
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raw (34) 


my/ (56) 


m z (36) 


m z (45) 


cross section (pb) 


o 


o 


o 


o 


6362 ± 0017 10~ 2 











1 


.1550 ±.0007 10~ 2 








1 





.5938 ±.0024 10~ 2 








1 


1 


.4144 ±.0038 10~ 3 





1 








.1441 ±.0002 10- 1 





1 





1 


.7142 ±.0045 10~ 3 





1 


1 





.8217 ±.0075 10~ 3 





1 


1 


1 


.2219 ±.0089 lO" 4 


1 











.1441 ±.0002 lO" 1 


1 








1 


.7194 ±.0045 10~ 3 


1 





1 





.8162 ±.0074 10~ 3 


1 





1 


1 


.2197 ±.0086 lO" 4 


1 


1 








.1318 ±.0003 10- 1 


1 


1 





1 


.6840 ±.0066 10^ 3 


1 


1 


1 





.1065 ±.0011 10- 3 


1 


1 


1 


1 


events 



Table 9: invariant mass analysis on e~ (3) u e (4) v e (5) e + (6) at \fs = 500 GeV. In 
the first four columns a number 1 is written when the invariant mass m v (ij) 
is inside the interval M v - 2T V < m v (ij) < My + 2T V (V= Z,W), other- 
wise the corresponding entry is 0. Additional cuts are E e + j e - > 20 GeV, 
| cos^ e+i e -| < 0.9, m(e+e-) > 30 GeV. 
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V~s = 


500 GeV 




V~s = 


1000 GeV 




V~s = 


2000 GeV 




(7(1) 


.2127 ± 


.0002 


io- 


-1 


.5885 ± 


.0010 


io- 


-2 


.2045 ± 


.0009 


io- 


-2 


a-(i) 


.1317 ± 


.0002 


10~ 


-1 


.3324 ± 


.0005 


10" 


-2 


.8707 ± 


.0022 


10" 


-3 


&tot (1) 


.1318 ± 


.0003 


10" 


-1 


.3339 ± 


.0022 


io- 


-2 


.8917 ± 


.0145 


io- 


-3 


cr(2) 


.6319 ± 


.0004 


io- 


-3 


.1548 ± 


.0002 


io- 


-3 


.4427 ± 


.0010 


10" 


-4 


a (2) 


.3921 ± 


.0003 


io- 


-3 


.9426 ± 


.0009 


io- 


-4 


.2589 ± 


.0003 


10" 


-4 


<5"tot (2) 


.4144 ± 


.0038 


10" 


-3 


.1009 ± 


.0015 


io- 


-3 


.2689 ± 


.0077 


10" 


-4 


(7(3) 


.1169 ± .0001 


io- 


-1 


.1201 ± 


.0001 


io- 


1 


.7250 ± 


.0018 


10" 


-2 


a (3) 


.8969 ± 


.0007 


io- 


-2 


.9483 ± 


.0009 


io- 


-2 


.5601 ± 


.0007 


io- 


-2 


0"tot (3) 


.1441 ± 


.0002 


io- 


-1 


.1161 ± 


.0004 


io- 


-1 


.6126 ± 


.0037 


io- 


-2 


cr (4) 


.1760 ± 


.0001 


io- 


-2 


.1159 ± 


.0001 


io- 


-2 


.6167 ± 


.0004 


io- 


-3 


a (4) 


.1379 ± 


.0001 


io- 


-2 


.9449 ± 


.0005 


io- 


-3 


.5103 ± 


.0004 


io- 


-3 


0"tot (4) 


.1550 ± 


.0007 


10" 


-2 


.1015 ± 


.0009 


io- 


-2 


.5281 ± 


.0060 


io- 


-3 


(7(5) 


.6840 ± 


.0013 


io- 


-2 


.1732 ± 


.0009 


io- 


1 


.3038 ± 


.0066 


10" 


-1 


6" (5) 


.5327 ± .0012 


io- 


-2 


.1382 ± 


.0007 


io- 


1 


.2233 ± 


.0031 


io- 


-1 


0"tot (5) 


.5938 ± 


.0024 


10" 


-2 


.1411 ± 


.0008 


io- 


-1 


.2275 ± 


.0031 


10" 


-1 




.4219 ± 


.0004 


10" 


-1 


.3653 ± 


.0010 


io- 


-1 


.4033 ± 


.0069 


io- 


-1 


Otot 


.6017 ± 


.0005 


10" 


-1 


.5099 ± 


.0011 


10" 


-1 


.4459 ± 


.0035 


10" 


-1 



Table 10: cross sections for e + e~v e v e in picobarns. a(i) are the cross sections 
for the processess (l)-(5) (see introduction) computed including only the 
signal diagrams, a (i) are the corresponding cross sections reconstructed using 
an invariant mass analysis (see caption of table |^). In o to t ah diagrams are 
included and & tot (i) are the results of the invariant mass analysis on a tot . 
E e+y e - > 20 GeV, | cos# e+i e -| < 0.9, m(e+e-) > 30 GeV. 
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^ = 500 GeV 


v^= 1000 GeV 


^ = 2000 GeV 


(7(1) 
(7(1) 

<5"tot (i) 


.4702 ±.0006 10" 1 
.2823 ±.0004 10" 1 
.2822 ±.0005 10" 1 


.1428 ±.0003 10" 1 
.7907 ±.0014 10~ 2 
.7920 ±.0020 10" 2 


.5460 ±.0025 10" 2 
.2272 ±.0006 10~ 2 
.2275 ±.0003 10" 2 


(7(3) 

a (3) 

&tot (3) 


.2961 ±.0002 10" 1 
.2279 ±.0002 10" 1 
.3516 ±.0004 10" 1 


.2953 ±.0003 10" 1 
.2348 ±.0002 10" 1 
.2836 ±.0003 10" 1 


.1817 ±.0001 10" 1 
.1417 ±.0001 10" 1 
.1559 ±.0001 10" 1 




.7663 ±.0008 10" 1 
.8273 ±.0007 10" 1 


.4381 ±.0006 10" 1 
.4520 ±.0004 10" 1 


.2363 ±.0004 10" 1 
.2192 ±.0001 10" 1 



Table 11: cross sections for e v e ud (pb). See caption of table [10 
E e - u j > 20 GeV, | cos^- u j| < 0.9, m(ud) > 30 GeV. 





^ = 500 GeV 


v^= 1000 GeV 


= 2000 GeV 


(7(2) 
(7(2) 
°tot (2) 


.7298 ±.0002 10" 3 
.4536 ± .0001 lO- 3 
.4765 ±.0008 10" 3 


.2074 ±.0001 10" 3 
.1301 ±.0001 10~ 3 
.1351 ±.0003 10" 3 


.6386 ±.0004 10" 4 
.3930 ±.0002 10~ 4 
.4015 ±.0012 10" 4 


(7(4) 

a (a) 

Otot (4) 


.1725 ±.0001 10" 2 
.1354 ±.0001 10~ 2 
.1361 ±.0001 10" 2 


.9679 ±.0002 10" 3 
.7874 ±.0002 10~ 3 
.7895 ±.0006 10" 3 


.4775 ±.0001 10" 3 
.3924 ± .0001 10" 3 
.3919 ±.0003 10" 3 


Otot 


.2455 ±.0001 10" 2 
.2833 ± .0001 10" 2 


.1175 ±.0001 10" 2 
.1398 ±.0001 10" 2 


.5414 ±.0001 10" 3 
.6471 ± .0004 10" 3 



Table 12: cross sections for e + e uu in picobarns. See caption of table [10 
E 



' (all particles) 



> 20 GeV, I cos 9 



(all particles) 



< 0.9, m(uu), m(e+e-) > 30 GeV. 
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^/i = 500 GeV 


y/s = 1000 GeV 


y/s = 2000 GeV 


(7(2) 

a (2) 

0"to* (2) 


.2165 ±.0001 10" 2 
.1334 ±.0001 10~ 2 
.1364 ±.0006 10" 2 


.5306 ±.0006 10" 3 
.3219 ±.0003 10~ 3 
.3273 ±.0045 10" 3 


.1507 ±.0003 10" 3 
.8848 ±.0013 10" 4 
.9275 ±.0282 10" 4 


(7(5) 
(7(5) 
^tot (5) 


.2242 ±.0005 10" 1 
.1749 ±.0004 lO" 1 
.1680 ±.0004 10" 1 


.5795 ±.0029 10" 1 
.4601 ±.0025 lO" 1 
.4554 ±.0025 10" 1 


.1331 ± .0324 
.7631 ±.0101 lO" 1 
.7781 ±.0104 10" 1 




.2459 ±.0005 10" 1 
.2374 ±.0004 10" 1 


.5848 ±.0029 10" 1 
.5663 ±.0027 10" 1 


.1333 ± .0324 
.9301 ±.0110 10" 1 



Table 13: Cross sections for v e v e uu (pb). See caption of table 
E u u > 20 GeV, | cos6 u a | < 0.9, m(uu) > 30 GeV. 
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